
Available online at www.sciencedirect.conn 



^^^^ • rv ■ Nuclear Physics B 

•# ScienceDirect proceedings 



Supplement 

ELSEVIER Nuclear Physics B Proceedings Supplement 00 (2013) l-[3] ^^^^^^^^= 



Supernova Detection in IceCube: Status and Future 

R. Bruijn for the IceCube Collaboration^ 

Laboratory for High Energy Physics, Ecole Polytechnique Federale, CH-1015, Lausanne, Switzerland 



Abstract 

The IceCube detector, located at the South Pole, is discussed as a detector for core collapse supernovae. The large 
flux of Ve from a Galactic supernova gives rise to Cherenkov light from positrons and electrons created in neutrino 
interactions which increase the overall count rate of the photomultipliers significantly. We will give an overview of 
the standard, count rate based, method for supernova detection and present the development of a novel technique. 
This technique uses coincident hits to extract additional information such as the average energy and spectral features. 
The potential of this technique increases with a higher sensor density, such as foreseen in projected extensions of 
IceCube/DeepCore . 

Keywords: Supernova, Neutrino, IceCube 



Introduction 

Many details of supernova explosions, like the exact 
mechanism, are still unknown, but it is clear that the 
largest part of the gravitational energy is released in 
the form of neutrinos with an energy in the order of 
10 MeV The detection of the neutrinos of a Galactic 
supernova can not only provide us with information on 
the dynamics of the explosion 1 1 1, but also on the nature 
of neutrinos such as the mass hierarchy. While pri- 
marily designed to detect astrophysical TeV neutrinos, 
the 5160 10" photomultipHers (PMTs) of the IceCube 
detector, instrumenting a cubic kilometer of ice, are 
used to detect MeV neutrinos from a Galactic super- 
nova |2, 3 1. The PMTs are housed in pressure-resistant 
glass spheres, together with read-out and digitization 
electronics, called DOMs. The DOMs are organized in 
vertical strings, reaching from 1450 to 2450 m depth. 
The 80 strings, each carrying 60 DOMs with a vertical 
spacing of 17.5 m, are arranged in a hexagonal grid 
with a leg lengths of about 125 m. Together with the 
surrounding strings, six additional strings carrying 
high quantum-efficiency PMTs, form the DeepCore 



sub-detector, which has a denser packing of about 60 m 
between strings and 7 m between DOMs on the same 
string. The primary channel for the detection of super- 
nova neutrinos is the inverse beta-decay interaction of 
Ve on protons. The produced positrons, which travel 
about 5 cm at an energy of 10 MeV, emit approximately 
2000 photons in the wavelength sensitivity range of 
the detector. The reference supernova model has a 8.8 
solar mass 0-Mg-Ne progenitor |4| and is considered 
conservative in terms of neutrino flux. This model is 
integrated over the first 3 seconds, optimizing for signal 
over background. We will first give an overview and 
status of supernova detection in IceCube and continue 
with a relatively new method which uses coincident hits 
to extract additional information such as average energy 
and other spectral features. This method is especially 
powerful with future projected dense detectors and 
results for a special geometry, called Deep and Dense 
Core (DDC) |5|, will be presented. This geometry 
consists of 24 strings with a 20 m spacing each equiped 
with 120 An sensitive DOMs separated by 3 m. 

Standard Method 

The detection of TeV neutrinos makes use of the 
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long distance space/time correlations between the 
photon hits caused by either a muon track or cascade 
initated by a neutrino. The low energy of supernova 
neutrinos in combination with the large separation of 
the DOMs make it impossible to reconstruct individual 
interactions. Instead, the detection technique relies on 
the vast amount of interactions in the instrumented 
volume, about 8 x 10^ for the reference model at 10 
kpc, which causes a significant increase in the count 
rate when considering all DOMs. Key to this is the 
low noise rate of the individual DOMs of about 540 
Hz, which is reduced to about 286 Hz after an artificial 
deadtime is applied after each hit, reducing the signal 
by 13%. The noise level is continuously monitored 
by a dedicated data-acquisition. It records scaler data 
in bins of 1.6384 ms. The data are later rebinned and 
synchronized in 2 ms bins. Following deployment, 
the noises rates show an exponential decrease which 
later levels out. A seasonal periodic behaviour with an 
amplitude of about 6% is seen in addition. The most 
likely collective rate deviation Aji from the baseline 
rate and its standard deviation (Ta// in a time interval 
are calculated in real-time with respect to two sliding 
windows of 300 s before and after the interval. Several 
values of At (0.5, 4 and 10 s) are used to take into 
account diff'erent models. Triggers are set at diff'erent 
levels of the significance ^ = Ajd/cr/^j^. When ^ exceeds 
6, an internal trigger is provided; when it exceeds 7.65, 
an alert is broadcasted on the Supernova Early Warning 
System (SNEWS). The significance values can be 
corrected for a correlation with the atmospheric muon 
rate. Significant triggers are caused by supemovae 
up to the Magellanic clouds (Fig. [T]). By considering 
the shape and amplitude of the lightcurves, assuming 
the models are known, diff'erent oscillation scenarios 
can be distinguished. Figure [2] shows the sensitivity 
to distinguish between normal and inverted hierarchy 
for diff'erent explosion models. The highest sensitivity 
is for scenarios in which a very massive star collapses 
into a quark star or black hole. In these scenarios, total 
enery and energies of the neutrinos are particularly 
large. 

Luminosity and energy 

A limitation of using only the singles rate comes from 
the entanglement of the luminosity and the energy 
distribution of the neutrinos |7|. The flux of neutrinos 
Oy^ at the detector from a supernova at distance d, can 
be described by 
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Figure 1: Expected significance versus distance assuming the 
Lawrence-Livermore model | 6 | for three oscillation scenarios. The 



density of the data points reflect the star distribution. 
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Figure 2: Number of standard deviations with which normal hierarchy 
and inverted hierarchy can be distinguished in at least 50% of all cases 

as function of supernova distance for several models. 

with L{t) describing the luminosity and fa(t),{E}(t)(Ev) the 
normaHzed energy distribution, which depends on the 
average energy {Ey)(t) and a 'pinch' factor a(t). For 
fa(t),{E)(t)(Ev), a commonly used parameterization is uti- 
lized: 
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The theoretical models provide L(t), {Ey}(t) and 
fa(t),{E}(t)(Ey) for each neutrino species. For a measured 
supernova, these are unknowns which have to be 
determined. In this work, we will assume that this 
formula can describe the neutrino energy distributions. 
From Eq[T] it is evident that an increase in luminosity 
can be compensated by a decrease in average energy. 
Furthermore, the detection probability depends on 
the neutrino energy; the cross-section and light yield 
increase quadratically and linearly with energy, re- 
spectively. To disentangle the luminosity, energy and 
other spectral features, additional energy sensitive 
observables are required. 
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Coincident Hits 

The probability of photons from a single interaction to 
reach different DOMs increases with the energy of the 
neutrino. We consider coincidences of 1 (single rate) 2 
or 3 DOMs, which can be nearest-neighbour, next- to- 
nearest-neighbour and on different strings. For 3 -fold 
coincidences 1, 2 or 3 strings can be involved. The 
different hit modes are sensitive to different parts of the 
energy spectrum. This sensitivity is determined by the 
number of DOMs and their geometry. The dependence 
of the probability is such that each extra DOM adds 
a factor of energy. This sensitivity of coincident hits 
to the neutrino energy has been recognized in fS] and 
is elaborated on in this work by including hit modes 
with more than two DOMs and a full Monte-Carlo 
simulation based on Geant4 with custom photon track- 
ing. The atmospheric muon background is taken into 
account, using a Gaisser parameterization of the flux. 
The uncorrelated background is greatly reduced by the 
use of a small coincidence gate of 150 ns for IceCube, 
100 ns for DeepCore and 50 for DDC. The correlated 
background from atmospheric muons can be reduced 
by a multiplicity cut. The required raw hit information 
is made available by an improved DAQ system which is 
currently being commissioned. Every hit is stored in a 
rotating filesystem, the hitspooling system. In case of a 
supernova alert, these data are transferred to permanent 
storage and are available for a detailed analysis, such as 
the coincident hit method. 

Determining energy and shape 

The energy dependent observables used are the ratios 
of the rates of the coincidence modes to the rate of 
single hits, i.e. the singles rate is used as an overall 
normalization, reducing systematic uncertainties. The 
ratio of the two-fold nearest-neighbour rate to the 
singles rate is indicated by rj^. Assuming the functional 
form of the energy distribution in Eq|2] with a fixed 
value Of, the energy dependence of r\^ on the average 
neutrino energy can be used to determine this energy 
from the rate |8|. The energy resolution obtained for 
IceCube and DDC using rj^, can be seen in Fig. [3] 
With IceCube, a resolution of 30% can be obtained 
for a supernova at 10 kpc. The power of the method, 
when used with a dense detector, is clear. With DDC, 
the resolution for the same conditions is improved to 
4%. By combining ratios, information on the spectral 
shape can be obtained. A first approach is constructing 
a test which includes ratios of all mentioned rates, 
and takes into account the correlations between the 
modes. This quantifies the difference between a 



measurement (simulation, in this work) and a chosen 
model (set of parameters a and {E)). A scan over the 
{{E), a) parameter space, presented in fig. |4j shows that 
the parameters can be constrained, however there is a 
degeneracv between them. 




Figure 3: Energy resolution, assuming a time integrated 8.8 solar 
mass O-Mg-Ne model with a = 2.84, Tor IceCube (solid) and DDC 

(dashed) at a distance of 5 (red), 10 (blue) and 20 (black) kpc. 




Figure 4: logO^^) contour as function of energy and of for a 8.8 O- 
Mg-Ne supernova at 10 kpc and the DDC detector geometry. True 

values indicated. 
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